Chapter 24 
Intermittency 



Sometimes They Come Back 
— Stephen King 

(R. Artuso, P. Dahlqvist, G. Tanner and P. Cvitanovic) 




N THE THEORY of chaotic dynamics developed so far we assumed that the evolu- 
tion operators have discrete spectra {zq,Zi,Z2, ■ ■ •) given by the zeros of 



k 

The assumption was based on the tacit premise that the dynamics is everywhere 
exponentially unstable. Real life is nothing like that - state spaces are gener- 
ically infinitely interwoven patterns of stable and unstable behaviors. The 
stable (in the case of Hamiltonian flows, integrable) orbits do not communicate 
with the ergodic components of the phase space, and can be treated by classical 
methods. In general, one is able to treat the dynamics near stable orbits as well 
as chaotic components of the phase space dynamics well within a periodic orbit 
approach. Problems occur at the borderline between chaos and regular dynamics 
where marginally stable orbits and manifolds present difficulties and still unre- 
solved challenges. 

We shall use the simplest example of such behavior - intermittency in 1- 
dimensional maps - to illustrate effects of marginal stability. The main message 
will be that spectra of evolution operators are no longer discrete, dynamical zeta 
functions exhibit branch cuts of the form 

i/^(z) = (--o + (i-zr(---), 

and correlations decay no longer exponentially, but as power laws. 
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24.1 Intermittency everywhere 

In many fluid dynamics experiments one observes transitions from regular behav- 
iors to behaviors where long time intervals of regular behavior ("laminar phases") 
are interrupted by fast irregular bursts. The closer the parameter is to the onset of 
such bursts, the longer are the intervals of regular behavior. The distributions of 
laminar- phase intervals are well described by power laws. 

This phenomenon is called intermittency, and it is a very general aspect of 
dynamics, a shadow cast by non-hyperbolic, marginally stable state space regions. 
Complete hyperbolicity assumed in (18.5) is the exception rather than the rule, 
and for almost any dynamical system of interest (dynamics in smooth potentials, 
billiards with smooth walls, the infinite horizon Lorentz gas, etc.) one encounters 
mixed state spaces with islands of stability coexisting with hyperbolic regions, 
see figure 24.1 andexample 7.7. Wherever stable islands are interspersed with 
chaotic regions, trajectories which come close to the stable islands can stay 'glued' 
for arbitrarily long times. These intervals of regular motion are interrupted by 
irregular bursts as the trajectory is re-injected into the chaotic part of the phase 
space. How the trajectories are precisely 'glued' to the marginally stable region is 
often hard to describe. What coarsely looks like a border of an island will under 
magnification dissolve into infinities of island chains of decreasing sizes, broken 
tori and bifurcating orbits, as illustrated in figure 24. 1 . 

Intermittency is due to the existence of fixed points and cycles of marginal 
stability (5.8), or (in studies of the onset of intermittency) to the proximity of a 
neaiiy mai^ginal complex or unstable orbits. In Hamiltonian systems intermittency 
goes hand in hand with the existence of (marginally stable) KAM tori. In more 
general settings, the existence of marginal or nearly marginal orbits is due to in- 
complete intersections of stable and unstable manifolds in a Smale horseshoe type 
dynamics (see figure 12.11). Following the stretching and folding of the invariant 
manifolds in time one will inevitably find state space points at which the stable and 
unstable manifolds are almost or exactly tangential to each other, implying non- 
exponential separation of nearby points in state space or, in other words, marginal 
stability. Under small parameter perturbations such neighborhoods undergo tan- 
gent bifurcations - a stable/unstable pair of periodic orbits is destroyed or created 
by coalescing into a marginal orbit, so the pruning which we shall encounter in 
chapter 12, and the intermittency discussed here are two sides of the same coin. section 12.4 
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Figure 24.2: A complete binary repeller with a 
marginal fixed point. 
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Figure 24.3: (a) A tent map trajectory, (b) A 
Farey map trajectory. 





How to deal with the full complexity of a typical Hamiltonian system with 
mixed phase space is a very difficult, still open problem. Nevertheless, it is pos- 
sible to learn quite a bit about intermittency by considering rather simple exam- 
ples. Here we shall restrict our considerations to 1-dimensional maps which in the 
neighborhood of a single marginally stable fixed point at x=0 take the form 



X ^ f{x) = x + O(x'+0 , 



(24.1) 



and are expanding everywhere else. Such a map may allow for escape, like the 
map shown in figure 24.2 or the dynamics may be bounded, hke the Farey map 



fix) 



x/{l-x) xe[0,l/2[ 
il-x)/x X€ [1/2,1] 



(24.2) 



Figure 24.3 compares a trajectory of the tent map (11.4) side by side with a 
trajectory of the Farey map. In a stark contrast to the uniformly chaotic trajectory 
of the tent map, the Farey map trajectory alternates intermittently between slow 
regular motion close to the marginally stable fixed point, and chaotic bursts. section 24.3.4 
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The presence of marginal stability has striking dynamical consequences: cor- 
relation decay may exhibit long range power law asymptotic behavior and diffu- 
sion processes can assume anomalous character. Escape from a repeller of the 
form figure 24.2 may be algebraic rather than exponential. In long time explo- 
rations of the dynamics intermittency manifests itself by enhancement of natural 
measure in the proximity of marginally stable cycles. 

The questions we shall address here are: how does marginal stability affect 
zeta functions or spectral determinants? And, can we deduce power law decays of 
correlations from cycle expansions? 

In example 23.5 we saw that marginal stability violates one of the conditions 
which ensure that the spectral determinant is an entire function. Already the sim- 
ple fact that the cycle weight 1/|1 - A^l in the trace (18.3) or the spectral determi- 
nant (19.3) diverges for marginal orbits with |Ap| = 1 tells us that we have to treat 
these orbits with care. 

In the following we will incorporate marginal stability orbits into cycle-expansions 
in a systematic manner. To get to know the difficulties lying ahead, we will start 
in sect. 24.2 with a piecewise linear map, with the asymptotics (24.1). We wiU 
construct a dynamical zeta function in the usual way without worrying too much 
about its justification and show that it has a branch cut singularity. We will cal- 
culate the rate of escape from our piecewise linear map and find that it is charac- 
terized by decay, rather than exponential decay, a power law. We will show that 
dynamical zeta functions in the presence of marginal stability can still be written 
in terms of periodic orbits, exactly as in chapters 17 and 22, with one exception: 
the mai^ginally stable orbits have to be explicitly excluded. This innocent looking 
step has far reaching consequences; it forces us to change the symbolic dynamics 
from a finite to an infinite alphabet, and entails a reorganization of the order of 
summations in cycle expansions, sect. 24.2.4. 

Branch cuts are typical also for smooth intermittent maps with isolated marginally 
stable fixed points and cycles. In sect. 24.3, we discuss the cycle expansions and 
curvature combinations for zeta functions of smooth maps tailored to intermit- 
tency. The knowledge of the type of singularity one encounters enables us to 
develop the efficient resummation method presented in sect. 24.3.1. 

Finally, in sect. 24.4, we discuss a probabilistic approach to intermittency that 
yields approximate dynamical zeta functions and provides valuable information 
about more complicated systems, such as billiards. 



24.2 Intermittency for pedestrians 

Intermittency does not only present us with a lai^ge repertoire of interesting dy- 
namics, it is also at the root of many son^ows such as slow convergence of cycle 
expansions. In order to get to know the kind of problems which arise when study- 
ing dynamical zeta functions in the presence of marginal stability we will consider 
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Figure 24.4: A piecewise linear intermittent map of 
(24.3) type: more specifically, the map piecewise lin- 
ear over intervals (24.9) of the toy example studied be- 
low, a = .5, b = .6, s = 1.0. 




an artfully concocted piecewise lineai^ model first. From there we will move on to 
the more general case of smooth intermittant maps, sect. 24.3. 



24.2.1 A toy map 

The Bernoulli shift map (23.6) is an idealized, but highly instructive, example 
of a hyperbolic map. To study intermittency we will now construct a likewise 
piecewise linear model, an intermittent map stripped down to its bare essentials. 

Consider a map x i-^ f{x) on the unit interval M = [0, 1] with two monotone 
branches 



fix) = 



foix) for X e Mo = [0, a] 
fi{x) for X € Ml = [b, 1] 



(24.3) 



The two branches are assumed complete, that is fo(Mo) = f\{Mi) = M. The map 
allows escape if a < b and is bounded if a = b (see figure 24.2 and figure 24.4). 
We take the right branch to be expanding and linear: 



1 



l-b 



(x-b). 



Next, we will construct the left branch in a way, which will allow us to model 
the intermittent behavior (24.1) near the origin. We chose a monotonically de- 
creasing sequence of points <7„ in [0,a] with q\ - a and <7„ — > as n ^ oo. 
This sequence defines a partition of the left interval Mq into an infinite number of 
connected intervals Mn, n>2 with 



Mn = ]qn,qn-\] and Alo ^ 1^ Mn- 



(24.4) 



n=2 



The map fo{x) is now specified by the following requirements 



foix) is continuous. 
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• fo{x) is linear on the intervals At„ for n > 2. 

• /o(^«) = qn~\, that is Mn = fo"*\[a, 1]) . 

This fixes the map for any given sequence {<7„). The last condition ensures the 
existence of a simple Markov partition. The slopes of the various linear segments 
are 

/oW = ^^ITrf^ = W forxeM„,n>3 

/oW - ^^^ITF^ = U\ fovxeM2 (24.5) 

with \Mn\ = <?n-i - Qn for n > 2. Note that we do not require as yet that the map 
exhibit intermittent behavior. 

We will see that the family of periodic orbits with code 10" plays a key role 
for intermittent maps of the form (24.1). An orbit 10" enters the intervals Mi — > 
Mn+\ — > Mn — > . . . ^ M2 successively and the family approaches the marginal 
stable fixed point at x = for n ^ 00. The stability of a cycle 10" for « > 1 is 
given by the chain rule (4.42), 

/ II 1 1 - a 

AlO" = /o(-^n+l)/o(-^«) ■ ■ • /o(-^2)/i(-Xl) ^ 7T7 [7— T ' (24-6) 

with X,- € Mi. 

The properties of the map (24.3) aie completely determined by the sequence 
{q„\. By choosing qn = 2"", for example, we recover the uniformly hyperbolic 
Bernoulli shift map (23.6). An intermittent map of the form (24.4) having the 
asymptotic behavior (24.1) can be constructed by choosing an algebraically de- 
caying sequence {qn} behaving asymptotically hke 

qn-4h^ (24.7) 

where s is the intermittency exponent in (24.1). Such a partition leads to intervals 
whose length decreases asymptotically like a power-law, that is, 

\Mn\ ~ ^ . (24.8) 

As can be seen from (24.6), the Floquet multipliers of periodic orbit families ap- 
proaching the marginal fixed point, such as the 10" family increase in turn only 
algebraically with the cycle length. 

It may now seem natural to construct an intermittent toy map in terms of a 
partition \M,i\ = l/n'"*"^^*, that is, a partition which follows (24.8) exactly. Such 
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a choice leads to a dynamical zeta function which can be written in terms of so- 
called Jonquiere functions (or polylogarithms) which arise naturally also in the 
context of the Farey map (24.2), and the anomalous diffusion of sect. 25.3. We remark 25.7 
will, however, not go along this route here; instead, we will engage in a bit of 
reverse engineering and construct a less obvious partition which will simplify the 
algebra considerably later without loosing any of the key features typical for inter- 
mittent systems. We fix the intermittent toy map by specifying the intervals A1„ 
in terms of Gamma functions according to 

r(n + m- \/s- 1) 

\Mn\=C^— '- for n>2, (24.9) 

i{n + m) 

where m = [l/s] denotes the integer part of l/s and C is a normalization constant 
fixed by the condition l^nl = - that is. 



C = a 



T{n-l/s) 



, nn + 1) 

n=m+l 



(24.10) 



Using Stirling's formula for the Gamma function 

r(z) ~ e-'-z''^'^ V2^(l + l/12z + ...), 

we verify that the intervals decay asymptotically like n~^^'*'^/''\ as required by the 
condition (24.8). 

Next, let us write down the dynamical zeta function of the toy map in terms 
of its periodic orbits, that is 



™-n(.-£) 



One may be tempted to expand the dynamical zeta function in terms of the binary 
symbolic dynamics of the map; we saw, however, in sect. 20.6 that such cycle ex- 
pansion converges extremely slowly. The shadowing mechanism between orbits 
and pseudo-orbits fails for orbits of the form 10" with stabilities given by (24.6), 
due to the marginal stability of the fixed point 0. It is therefore advantageous to 
choose as the fundamental cycles the family of orbits with code 10" or, equiva- 
lently, switch from the finite (binary) alphabet to an infinite alphabet given by 



10"-i ^ n. 



Due to the piecewise-linear form of the map which maps intervals Mn exactly 
onto Mn-i, all periodic orbits entering the left branch at least twice ai^e canceled 
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exactly by pseudo cycles, and the cycle expanded dynamical zeta function depends 
only on the fundamental series 1, 10, 100, . . .: 



The fundamental term (20.7) consists here of an infinite sum over algebraically 
decaying cycle weights. The sum is divergent for \z\ > 1. We will see that this 
behavior is due to a branch cut of 1 starting at z = 1 . We need to find analytic 
continuations of sums over algebraically decreasing terms in (24.11). Note also 
that we omitted the fixed point in the above Euler product; we will discussed 
this point as well as a proper derivation of the zeta function in more detail in 
sect. 24.2.4. 



24.2.2 Branch cuts 

Starting from the dynamical zeta function (24.11), we first have to worry about 
finding an analytical continuation of the sum for |z| > 1 . We do, however, get this 
part for free here due to the particular choice of interval lengths made in (24.9). 
The sum over ratios of Gamma functions in (24.11) can be evaluated analytically 
by using the following identities valid for I/5 = a > (the famed binomial 
theorem in disguise), 

• a non-integer 





(24.11) 




(24.12) 



• a integer 



a 




(24.13) 



+ 




with 
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In order to simplify the notation, we restrict the intermittency parameter to the 
range 1 < l/s <2 with [l/s] = m = 1. All what follows can easily be generalized 
to arbitrary s > using equations (24.12) and (24.13). The infinite sum in (24.11) 
can now be evaluated with the help of (24.12) or (24.13), that is, 

y Tin- l/s) „ ^ r r(-i) [(1 - zy'' - 1 + iz] for 1 < I/5 < 2; 
^r{n+l)^ \ (1 -z)log(l -z) + z ■ for s^\. 

The normalization constant C in (24.9) can be evaluated explicitly using (24.10) 
and the dynamical zeta function can be given in closed form. We obtain for 1 < 

\/s < 2 

\/az) ^\-{l-b)z- ,/ ■ ((1 - z)"' - 1 + -z] . (24.14) 
l/s - I I - a\ s ) 

and for s = I, 

l/az) = l-(,l-b)z- ^ ((1 - z)log(l -z) + z). (24.15) 

It now becomes clear why the particular choice of intervals Mn made in the last 
section is useful; by summing over the infinite family of periodic orbits 0" 1 ex- 
plicitly, we have found the desired analytical continuation for the dynamical zeta 
function for |z| > 1 . The function has a branch cut starting at the branch point z = 1 
and running along the positive real axis. That means, the dynamical zeta function 
takes on different values when approaching the positive real axis for Re z > 1 from 
above and below. The dynamical zeta function for general > takes on the form 

l/^(z) - 1 - (1 - - ((1 - z)"' - gAz)) (24.16) 

gsW l-az'" ^ ^ ' 

for non-integer s with m - [I/5] and 

l/^(z) - 1 - (1 - z)'"log(l - z) - g„,(z))(24.I7) 

for \js = m integer and gsiz) are polynomials of order m = [l/s] which can 
be deduced from (24.12) or (24.13). We thus find algebraic branch cuts for non 
integer intermittency exponents l/s and logarithmic branch cuts for l/s integer. 
We will see in sect. 24.3 that branch cuts of that form are generic for 1-dimensional 
intermittent maps. 

Branch cuts are the all important new feature of dynamical zeta functions due 
to intermittency. So, how do we calculate averages or escape rates of the dynamics 
of the map from a dynamical zeta function with branch cuts? We take 'a learning 
by doing' approach and calculate the escape from our toy map for a < b. 



inter - 19iiov2012 



ChaosBook.org versioiiI4.4.1, Apr 21 2013 



CHAPTER 24. INTERMITTENCY 



488 



Figure 24.5: The survival probability r„ calcu- 
lated by contour integration; integrating (24.18) 
inside the domain of convergence \z\ < 1 (shaded 
area) of l/i^iz) in periodic orbit representation 
yields (18.26). A deformation of the contour 
(dashed line) to a larger circle gives contribu- 
tions from the poles and zeros (x) of l/((z) be- 
tween the two circles. These are the only contribu- 
tions for hyperbolic maps (a), for intermittent sys- 
tems additional contributions arise, given by the 
contour 7„„ running along the branch cut (b). 



(a) 




24.2.3 Escape rate 

Our starting point for the calculation of the fraction of survivors after n time steps, 
is the integral representation (19.19) 



where the contour encircles the origin in the clockwise direction. If the contour 
lies inside the unit circle \z\ = 1, we may expand the logarithmic derivative of 
(~^iz) as a convergent sum over all periodic orbits. Integrals and sums can be 
interchanged, the integrals can be solved term by term, and the formula (18.26) 
is recovered. For hyperbolic maps, cycle expansion methods or other techniques 
may provide an analytic extension of the dynamical zeta function beyond the lead- 
ing zero; we may therefore deform the original contour into a larger circle with ra- 
dius R which encircles both poles and zeros of ^~^{z), see figure 24.5 (a). Residue 
calculus turns this into a sum over the zeros Za and poles z^ of the dynamical zeta 
function, that is 



where the last term gives a contribution from a large circle y^. We thus find 
exponential decay of r„ dominated by the leading zero or pole of ^"^(z). 

Things change considerably in the intermittent case. The point z = 1 is a 
branch cut singularity and there exists no Taylor series expansion of ("^ ai^ound 
z = 1. Second, the path deformation that led us to (24.19) requires more care, as it 
must not cross the branch cut. When expanding the contour to large \z\ values, we 
have to deform it along the branch Re (z) > 1 , Im (z) = encircling the branch cut 
in anti-clockwise direction, see figure 24.5 (b). We will denote the detour around 
the cut as ycut- may write symbolically 




(24.18) 




(24.19) 




zeros poles 
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where the sums include only the zeros and the poles in the area enclosed by the 
contours. The asymptotics is controlled by the zero, pole or cut closest to the 
origin. 

Let us now go back to our intermittent toy map. The asymptotics of the sur- 
vival probability of the map is here governed by the behavior of the integrand 
^ log in (24.18) at the branch point z = 1. We restrict ourselves again to the 
case I < 1/ s < 2 first and write the dynamical zeta function (24.14) in the form 

l/^(z) = ao + ai(\-z) + bo(\ - z)^'' = G(l - z) 

and 

b - a , a I — b 



ao 



I — a I — 1/ s I — a 

Setting M = 1 - z, we need to evaluate 



(fi ( 1 - u)-" 4- log G{u)du (24.20) 

2m Jy^„, du 



where -ydt goes around the cut (i.e., the negative u axis). Expanding the integrand 

d_ 

du 



4- \ogG{u) = G'{u)/G{u) in powers of u and w'^* at m = 0, one obtains 



^logG(«) = ^ + -^u''-^-' + 0{u). (24.21) 
du ao s ao 



The integrals along the cut may be evaluated using the general formula 



— £ u"{l- u)~"du = — ]1 L(i + 0(1 /n)) (24.22) 



which can be obtained by deforming the contour back to a loop around the point 
u = I, now in positive (anti-clockwise) direction. The contour integral then picks 
up the («-l)st term in the Taylor expansion of the function u" atu = 1, cf. (24. 12). 
For the continuous time case the corresponding formula is 



Plugging (24.21) into (24.20) and using (24.22) we get the asymptotic result 



bal 1 1 a l-b 1 1 

r„ ~ — = — . (24.24) 

ao sY{l - l/s)n^/'^ s - I b - aY{\ - \/s) n^'' 
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Figure 24.6: The asymptotic escape from an inteiTnit- 
tent repeller is a power law. Normally it is preceded 
by an exponential, which can be related to zeros close 
to the cut but beyond the branch point z = 1, as in 
figure 24.5 (b). 




200 400 600 800 1000 



We see that, asymptotically, the escape from an intermittent repeller is described 
by power law decay rather than the exponential decay we are familiar with for 
hyperbolic maps; a numerical simulation of the power-law escape from an inter- 
mittent repeller is shown in figure 24.6. 

For general non-integer l/s > 0, we write 
l/^(z) = A{u) + {u)^l'B{u) = G{u) 

with u = \ - z and A{u), B{u) are functions analytic in a disc of radius 1 around 
u -0. The leading terms in the Taylor series expansions of A{u) and B{u) are 

b — a a \ —b 

flo = -; , bo = 



1-a' g,(l)l-fl'' 



see (24.16). Expanding ^logG(M) around u = 0, one again obtains leading or- 
der contributions according to (24.21) and the general result follows immediately 
using (24.22), that is, 

a l-b 1 1 

^" ~ sg,{l)b-ar{l - l/s)n^/' ' 

Applying the same arguments for integer intermittency exponents l/s = m, one 
obtains 

(24.26) 

sgmW b-an'" 

So far, we have considered the survival probability for a repeller, that is we 
assumed a < b. The formulas (24.25) and (24.26) do obviously not apply for the 
case a = b, that is, for the bounded map. The coefficient aQ = (b - a)/{l - a) 
in the series representation of G{u) is zero, and the expansion of the logarithmic 
derivative of G(u) (24.21) is no longer valid. We get instead 



— logG(M) - <^ "), , < 
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assuming non-integer 1 for convenience. One obtains for the survival probabil- 
ity. 




1 + Oin^-^'"') s < 1 
l/s + 0{n^''~^) s>l 



For s > \, this is what we expect. There is no escape, so the survival probability 
is equal to 1, which we get as an asymptotic result here. The result for 5 > 1 is 
somewhat more worrying. It says that F,, defined as sum over the instabilities of 
the periodic orbits as in (22.18) does not tend to unity for large n. However, the 
case 5 > 1 is in many senses anomalous. For instance, the invariant density cannot 
be normalized. It is therefore not reasonable to expect that periodic orbit theories 
will work without complications. 

24.2.4 Why does it work (anyway)? 

Due to the piecewise linear nature of the map constructed in the previous section, 
we had the nice property that interval lengths did exactly coincide with the inverse 
of the stability of periodic orbits of the system, that is 



There is thus no problem in replacing the survival probability F„ given by (1.2), 
(22.2), that is the fraction of state space M surviving n iterations of the map. 



by a sum over periodic orbits of the form (18.26). The only orbit to worry about is 
the mai^ginal fixed point itself which we excluded from the zeta function (24. 11). 

For smooth intermittent maps, things are less clear and the fact that we had to 
prune the marginal fixed point is a warning sign that interval estimates by periodic 
orbit stabilities might go horribly wrong. The derivation of the survival probability 
in terms of cycle stabilities in chapter 22 did indeed rely heavily on a hyperbolicity 
assumption which is clearly not fulfilled for intermittent maps. We therefore have 
to carefully reconsider this derivation in order to show that periodic orbit formulas 
are actually valid for intermittent systems in the first place. 

We will for simplicity consider maps, which have a finite number of say s 
branches defined on intervals Ms and we assume that the map maps each inter- 
val Ms onto M, that is f{Ms) = M. This ensures the existence of a complete 
symbolic dynamics - just to make things easy (see figure 24.2). 



\Mn\ - l/|Aior 




(«) 
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The generating partition is composed of the domains Ms. The ?ith level parti- 
tion C^"^ = {Mi] can be constructed iteratively. Here /'s are words / = S2S2 ■ ■ .s„ 
of length n, and the intervals Mi are constructed recursively 

Msj^f;\Mj), (24.27) 

where sj is the concatenation of letter s with word j of length nj < n. 

In what follows we will concentrate on the survival probability r„ , postponing 
other quantities of interest, such as averages, to later considerations. In establish- 
ing the equivalence of the survival probability and the periodic orbit formula for 
the escape rate for hyperbolic systems we have assumed that the map is expand- 
ing, with a minimal expansion rate \f'{x)\ > Amin > 1- This enabled us to bound 
the size of every survivor strip Mi by (22.6), the stability A,- of the periodic orbit / 
within the Mi, and bound the survival probability by the periodic orbit sum (22.7). 



The bound (22.6) 

'|A,| |M| -|A,| 

relies on hyperbolicity, and is thus indeed violated for intermittent systems. The 
problem is that now there is no lower bound on the expansion rate, the minimal 
expansion rate is A^in = 1 . The survivor strip Mo" which includes the marginal 
fixed point is thus completely overestimated by l/|Ao"| = 1 which is constant for 
all n. exercise 19.7 

However, bounding survival probability strip by strip is not what is required 
for establishing the bound (22.7). For intermittent systems a somewhat weaker 
bound can be established, saying that the average size of intervals along a periodic 
orbit can be bounded close to the stability of the periodic orbit for all but the 
interval M(y. The weaker bound applies to averaging over each prime cycle p 
separately 

C,J-<lyk^<C2^, (24.28) 

where the word / represents a code of the periodic orbit p and all its cyclic permu- 
tations. It can be shown that one can find positive constants Ci, C2 independent 
of p. Summing over all periodic orbits leads then again to (22.7). 

To study averages of multiplicative weights we follow sect. 17. 1 and introduce 
a state space observable a{x) and the integrated quantity 

n-l 



A"(x)-2]a(/(x)). 



k=0 
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o 

Figure 24.7: Transition graph corresponding to the al- 
phabet {0*'-' 1;(), k>l] 




This leads us to introduce the generating function (17.10) 

where (.) denote some averaging over the distribution of initial points, which we 
choose to be uniform (rather than the a priori unknown invariant density). Again, 
all we have to show is, that constants Ci, C2 exist, such that 



C,^<l2]^f c.^, (24.29) 

|Ap| rip \M\ Jmq |Ap| 



is valid for all p. After performing the above average one gets 

<^ir„(yS)<^ r < C2r„0S), (24.30) 



with 



and a dynamical zeta function can be derived. In the intermittent case one can 
expect that the bound (24.29) holds using an averaging argument similar to the 
one discussed in (24.28). This justifies the use of dynamical zeta functions for 
intermittent systems. 

One lesson we should have learned so far is that the natural alphabet to use 
is not {0, 1} but rather the infinite alphabet {0*^~'l,0; k > 1}. The symbol 
occurs unaccompanied by any I's only in the marginal fixed point which is 
disconnected from the rest of the transition graph, see figure 24.7. chapter 12 

What happens if we remove a single prime cycle from a dynamical zeta func- 
tion? In the hyperbolic case such a removal introduces a pole in the 1 and slows 
down the convergence of cycle expansions. The heuristic interpretation of such a 
pole is that for a subshift of finite type removal of a single prime cycle leads to 
unbalancing of cancellations within the infinity of of shadowing pairs. Neverthe- 
less, removal of a single prime cycle is an exponentially small perturbation of the 
trace sums, and the asymptotics of the associated trace formulas is unaffected. chapter 23 
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In the intermittent case, the fixed point does not provide any shadowing , 
and a statement such as 



is meaningless. It seems therefore sensible to take out the factor (l - to) = I - z 
from the product representation of the dynamical zeta function (19.15), that is, to 
consider a pruned dynamical zeta function ll^interiz) defined by 



We saw in the last sections, that the zeta function ll^interiz) has all the nice prop- 
erties we know from the hyperbolic case, that is, we can find a cycle expansion 
with - in the toy model case - vanishing curvature contributions and we can calcu- 
late dynamical properties like escape after having understood, how to handle the 
branch cut. But you might still be wonied about leaving out the extra factor 1 - z 
all together. It turns out, that this is not only a matter of convenience, omitting 
the marginal cycle is a dire necessity. The cycle weight = 1 overestimates 
the corresponding interval length of Mq" in the partition of the state space M by 
an increasing amount thus leading to wrong results when calculating escape. By 
leaving out the cycle (and thus also the Alo" contribution), we are guaranteed to 
get at least the right asymptotical behavior. 

Note also, that if we are working with the spectral determinant (19.3), given 
in product form as 



for intermittent maps the marginal stable cycle has to be excluded. It introduces 
an (unphysical) essential singularity at z = 1 due the presence of a factor (1 - z)°° 
stemming from the cycle. 



24.3 Intermittency for cyclists 

Admittedly, the toy map is what is says - a toy model. The piece wise linear- 
ity of the map led to exact cancellations of the curvature contributions leaving 
only the fundamental terms. There are still infinitely many orbits included in the 
fundamental term, but the cycle weights were chosen in such a way that the zeta 
function could be written in closed form. For a smooth intermittent map this all 
will not be the case in general; still, we will argue that we have already seen al- 
most all the fundamentally new features due to intermittency. What remains are 
technicalities - not necessarily easy to handle, but nothing very surprise any more. 



Ai.qa+i ~ Ai.QtAo, 



\iaz) = {l-z)\IUer{z). 
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In the following we will sketch, how to make cycle expansion techniques work 
for general 1 -dimensional maps with a single isolated marginal fixed point. To 
keep the notation simple, we will consider two-branch maps with a complete bi- 
nary symbolic dynamics as for example the Farey map, figure 24.3, or the repeller 
depicted in figure 24.2. We again assume that the behavior near the fixed point is 
given by (24.1). This implies that the stability of a family of periodic orbits ap- 
proaching the marginally stable orbit, as for example the family 10", will increase 
only algebraically, that is we find again for large n 

1 1 

Aio" ~ n^^^/' ' 

where s denotes the intermittency exponent. 

When considering zeta functions or trace formulas, we again have to take out 
the marginal orbit 0; periodic orbit contributions of the form ?o"i are now unbal- 
anced and we arrive at a cycle expansion in terms of infinitely many fundamental 
terms as for our toy map. This corresponds to moving from our binary symbolic 
dynamics to an infinite symbolic dynamics by making the identification 

10"-' ^ n; 10"-' 10'"-' ^ nm; 10"-' 10"'-' 10*=-' ^ nmk; . . . 



see also table 24. 1. The topological length of the orbit is thus no longer determined 
by the iterations of our two-branch map, but by the number of times the cycle 
goes from the right to the left branch. Equivalently, one may define a new map, 
for which all the iterations on the left branch are done in one step. Such a map is 
called an induced map and the topological length of orbits in the infinite alphabet 
corresponds to the iterations of this induced map. exercise 12.1 

For generic intermittent maps, curvature contributions in the cycle expanded 
zeta function will not vanish exactly. The most natural way to organize the cycle 
expansion is to collect orbits and pseudo orbits of the same topological length 
with respect to the infinite alphabet. Denoting cycle weights in the new alphabet 
as tnm... = fio"-'iO'"-i...' oris obtains 



- Y\{l-tp)^l-Y,ce (24.32) 

CO CO CO ^ 

~ 1 ~ ^ ] tfi ~ ^ ] ^ ] "^(.fmn ~ tf,jt„) 
n=i m=l n=l 

CO CO CO . . ^ DO DO DO CO 

k=l m=\ n=\ l=\ k=\ m=l n=l 



The first sum is the fundamental term, which we have already seen in the toy 
model, (24.11). The curvature terms c„ in the expansion are now e-fold infinite 
sums where the prefactors take care of double counting of prime periodic orbits. 
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Table 24.1: Infinite alphabet versus the original binary alphabet for the shortest periodic 
orbit families. Repetitions of prime cycles (11 = 1^,0101 = 01^, . . .) and their cyclic 
repeats (110 = 101,1110 = 1101,...) are accounted for by cancelations and combination 
factors in the cycle expansion (24.32). 



oo - alphabet 






binary alphabet 








n = 1 


n = 2 


« = 3 


n = 4 


« = 5 


1 -cycles 


n 


1 


10 


100 


1000 


10000 


2-cycles 


mn 














In 


11 


110 


1100 


11000 


110000 




2n 


101 


0101 


10100 


101000 


1010000 




3n 


1001 


10010 


100100 


1001000 


10010000 




An 


10001 


100010 


1000100 


10001000 


100010000 


3 -cycles 


kmn 














\\n 


111 


1110 


11100 


111000 


1110000 




I2n 


1101 


11010 


110100 


1101000 


11010000 




13n 


11001 


110010 


1100100 


11001000 


110010000 




21n 


1011 


10110 


101100 


1011000 


10110000 




22n 


10101 


101010 


1010100 


10101000 


101010000 




23n 


101001 


1010010 


10100100 


101001000 


1010010000 




31n 


10011 


100110 


1001100 


10011000 


100110000 




32« 


100101 


1001010 


10010100 


100101000 


1001010000 




33« 


1001001 


10010010 


100100100 


1001001000 


10010010000 



Let us consider the fundamental term first. For generic intermittent maps, we 
can not expect to obtain an analytic expression for the infinite sum of the form 

oo 

m = Y,hnZ". (24.33) 

n=0 

with algebraically decreasing coefficients 

hn with a > 

n" 

To evaluate the sum, we face the same problem as for our toy map: the power 
series diverges for z > 1, that is, exactly in the 'interesting' region where poles, 
zeros or branch cuts of the zeta function are to be expected. By carefully subtract- 
ing the asymptotic behavior with the help of (24. 12) or (24. 13), one can in general 
construct an analytic continuation of /(z) around z = 1 of the form 



/(z) ~ A{z) + (I - zf'^ B{z) atn (24.34) 

/(z) ~ A(z) + (1 -z)''~'ln(l -z) a en, 

where A(z) and B(z) are functions analytic in a disc around z = 1. We thus again 
find that the zeta function (24.32) has a branch cut along the real axis Rez > 1- 
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From here on we can switch to auto-pilot and derive algebraic escape, decay of 
correlation and all the rest. We find in particular that the asymptotic behavior 
derived in (24.25) and (24.26) is a general result, that is, the survival probability 
is given asymptotically by 

r„ ~ C^TT (24.35) 

for all 1-dimensional maps of the form (24.1). We have to work a bit harder if 
we want more detailed information like the prefactor C, exponential precursors 
given by zeros or poles of the dynamical zeta function or higher order corrections. 
This information is buried in the functions A{z) and B{z) or more generally in the 
analytically continued zeta function. To get this analytic continuation, one may 
follow either of the two different strategies which we will sketch next. 



24.3.1 Resummation 



One way to get information about the zeta function near the branch cut is to de- 
rive the leading coefficients in the Taylor series of the functions A(z) and B{z) in 
(24.34) at z = 1. This can be done in principle, if the coefficients /i„ in sums like 
(24.33) are known (as for our toy model). One then considers a resummation of 
the form 

oo oo oo 

2] V = 2 - zy + (1 - zf-' ^/l - 2)'' (24.36) 
y=o i=o j=o 

and the coefficients aj and bj are obtained in terms of the hj's by expanding (1 -zY 
and (1 - z)^^''"' on the right hand side around z = using (24.12) and equating 
the coefficients. 

In practical calculations one often has only a finite number of coefficients 
hj, < j < N, which may have been obtained by finding periodic orbits and 
their stabilities numerically. One can still design a resummation scheme for the 
computation of the coefficients aj and bj in (24.36). We replace the infinite sums 
in (24.36) by finite sums of increasing degrees and nt, and require that 

'la It, N 

J] aid - zi + (1 - zf-' - = Z ^^"^ + ^^"-"^'^ ■ (^'^•^"^^ 

(=0 (=0 j=0 

One proceeds again by expanding the right hand side around z - 0, skipping all 
powers z^^^ and higher, and then equating coefficients. It is natural to require that 
\nb + a - \ - na\ < 1, so that the maximal powers of the two sums in (24.37) are 
adjacent. If one chooses ria + rih + 2 = N + I, then, for each cutoff length A'^, the 
integers and Uh are uniquely determined from a linear system of equations. The 
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price we pay is that the so obtained coefRcients depend on the cutoff A'^. One can 
now study convergence of the coefRcients aj, and bj, with respect to increasing 
values of N, or various quantities derived from aj and bj. Note that the leading 
coefficients and bo determine the prefactor C in (24.35), cf. (24.24). The re- 
summed expression can also be used to compute zeros, inside or outside the radius 
of convergence of the cycle expansion ^ hjz^. 

The scheme outlined in this section tacitly assumes that a representation of 
form (24.34) holds in a disc of radius 1 around z = 1. Convergence is improved 
further if additional information about the asymptotics of sums like (24.33) is used 
to improve the ansatz (24.36). 



24.3.2 Analytical continuation by integral transformations 

We will now introduce a method which provides an analytic continuation of sums 
of the form (24.33) without explicitly relying on an ansatz (24.36). The main 
idea is to rewrite the sum (24.33) as a sum over integrals with the help of the 
Poisson summation formula and find an analytic continuation of each integral by 
contour deformation. In order to do so, we need to know the n dependence of 
the coefficients /i„ = h(n) explicitly for all n. If the coefficients are not known 
analytically, one may proceed by approximating the large n behavior in the form 

h(n) - n'^iCi + C2n~^ + ...), ntO, 

and determine the constants C, numerically from periodic orbit data. By using the 
Poisson resummation identity 

oo oo 

^ S(x-n)^ ^ exp(2nimx) , (24.38) 

n=-oo m=-oo 

we may write the sum as (24.33) 
1 . r*°^ 

f{z) = -hm+Y, I dxe^"'"'''h{x)z\ (24.39) 

m=—oo ^ 

The continuous variable x con^esponds to the discrete summation index n and it 
is convenient to write z = r exp(/cr) from now on. The integrals are still not con- 
vergent for r > 0, but an analytical continuation can be found by considering the 
contour integral, where the contour goes out along the real axis, makes a quarter 
circle to either the positive or negative imaginary axis and goes back to zero. By 
letting the radius of the circle go to infinity, we essentially rotate the line of inte- 
gration from the real onto the imaginary axis. For the m = term in (24.39), we 
transform x ix and the integral takes on the form 

dxh{x)r'e''"^ = i \ dxh{ix)r'-'e-"^. 
Jo Jo 
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The integrand is now exponentially decreasing for all r > and cr or In. The 
last condition reminds us again of the existence of a branch cut at Rez > 1. By 
the same technique, we find the analytic continuation for all the other integrals in 
(24.39). The real axis is then rotated according to — > sign(m)/;c where sign(m) 
refers to the sign of m. 



■\m\±<T) 



Jo Jo 



Changing summation and integration, we can carry out the sum over \m\ explicitly 
and one finally obtains the compact expression 



f(z) = -/<0) + / r dxhiix)r'-'e-"^ (24.40) 
2 Jo 

Xoo -2;r.v 
dx J— ^ [Kixy-'e-"^ - K-ix)r-'''e'"^] . 



The transformation from the original sum to the two integrals in (24.40) is exact 
for r < 1 , and provides an analytic continuation for r > 0. The expression (24.40) 
is especially useful for an efficient numerical calculations of a dynamical zeta 
function for |z| > 1, which is essential when searching for its zeros and poles. 



24.3.3 Curvature contributions 



So far, we have discussed only the fundamental term in (24.32), and 

showed how to deal with such power series with algebraically decreasing coef- 
ficients. The fundamental term determines the main structure of the zeta function 
in terms of the leading order branch cut. Corrections to both the zeros and poles 
of the dynamical zeta function as well as the leading and subleading order terms 
in expansions like (24.34) ai^e contained in the curvature terms in (24.32). The 
first curvature correction is the 2-cycle sum 



DO OO ^ 

m=l n=\ 



tmtn) ^ 



with algebraically decaying coefficients which again diverge for \z\ > 1. The 
analytically continued curvature terms have as usual branch cuts along the positive 
real z axis. Our ability to calculate the higher order curvature terms depends on 
how much we know about the cycle weights ?„„,. The form of the cycle stability 
(24.6) suggests that t,„„ decrease asymptotically as 



(nm) 



l + l/s 



(24.41) 
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for 2-cycles, and in general for ^-cycles as 



1 



t, 



{m\m2 ■ ■ ■ m„)'+i''' 



hi 



If we happen to know the cycle weights tmim2...m„ analytically, we may proceed as 
in sect. 24.3.2, transform the multiple sums into multiple integrals and rotate the 
integration contours. 

We have reached the edge of what has been accomplished so far in computing 
and what is worth the dynamical zeta functions from periodic orbit data. In the 
next section, we describe a probabilistic method appUcable to intermittent maps 
which does not rely on periodic orbits. 

24.3.4 Stability ordering for intermittent flows 



Longer but less unstable cycles can give larger contributions to a cycle 
expansion than short but highly unstable cycles. In such situations, truncation by 
length may require an exponentially large number of very unstable cycles before 
a significant longer cycle is first included in the expansion. This situation is best 
illustrated by intermittent maps. The simplest of these is the Farey map 



For the Farey map, the symbolic dynamics is of complete binary type, so the 
lack of shadowing is not due to the lack of a finite grammar, but rather to the 
intermittency caused by the existence of the marginal fixed point xq = 0, for which 
the stability multiplier is Aq = 1 . This fixed point does not participate directly in 
the dynamics and is omitted from cycle expansions. Its presence is, however, 
very much felt instead in the stabilities of neighboring cycles with n consecutive 
iterates of the symbol 0, whose stability falls of only as A ~ n^, in contrast to 
the most unstable cycles with n consecutive I's, which are exponentially unstable. 



The symbolic dynamics is of complete binary type. A quick count in the style 
of sect. 15.7.2 leads to a total of 74,248,450 prime cycles of length 30 or less, 
not including the marginal point xq = 0. Evaluating a cycle expansion to this 
order is an impressive computational feat. However, stability of the least unstable 
cycle omitted is roughly Ajoso ~ 30^ = 900, so it yields a 0.1% correction. The 
situation may be much worse than this estimate suggests, because the next 10^' 
cycle contributes a similar amount, and could easily reinforce the error. Adding 
up all such omitted terms, we arrive at an estimated error of about 3%, for a cycle- 
length truncated cycle expansion based on more than 10^ pseudo-cycle terms! On 
the other hand, if one truncates by stability at Amax - 3000, only 409 prime cycles 
suffice to attain the same accuracy of about a 3% error, figure 24.8. 





< X < 1/2 
1/2 < x< 1 



(24.42) 



|Aor|~[(V5 + l)/2] 
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Figure 24.8: Comparison of cycle expansion tran- 
cation schemes for the Farey map (24.42); the 
deviation of the truncated cycles expansion for 
|l/^iv(0)| from the exact flow conservation value 
l/<r(0) = is a measure of the accuracy of the 
truncation. The jagged line is the logarithm of 
the stability ordering truncation error; the smooth 
line is smoothed according to sect. 20.6.2; the di- 
amonds indicate the error due to the topological 
length truncation, with the maximal cycle length N 
shown. They are placed along the stability cutofl" 
axis at points determined by the condition that the 
total number of cycles is the same for both trunca- 
tion schemes. 
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As the Farey map maps the unit interval onto itself, the leading eigenvalue 
of the Penon-Frobenius operator should equal = 0, so 1/^(0) = 0. Devia- 
tions from this exact result give an indication of the convergence of a given cycle 
expansion. Errors corresponding to different truncation schemes are indicated in 
figure 24.8. We see that topological length truncation schemes are hopelessly bad 
in this case; stability length truncations are somewhat better, but still rather bad. 
In simple cases like this one, where intermittency is caused by a single marginal 
fixed point, convergence can be improved by going to infinite alphabets. 



^^Or So far we have focused on 1-d models as the simplest setting in which to 
investigate dynamical implications of marginal fixed points. We now take an al- 
together different track and describe how probabilistic methods may be employed 
in order to write down approximate dynamical zeta functions for intermittent sys- 
tems. 

We will discuss the method in a very general setting, for a flow in arbitrary 
dimension. The key idea is to introduce a surface of section P such that all tra- 
jectories traversing this section will have spent some time both neai" the marginal 
stable fixed point and in the chaotic phase. An important quantity in what follows 
is (3.5), the first return time t{x), or the time of flight of a trajectory starting in 
X to the next return to the surface of section "P. The period of a periodic orbit p 
intersecting the P section rip times is 



where /(x) is the Poincare map, and Xp e ;P is a periodic point. The dynamical 



24.4 BER zeta functions 





1=0 
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zeta function (19.15) 



l/az, s,p) = II 1 — , - 2] «(/ (^p))' (24-43) 

P ^ ' p' ' k=0 



chapter 17 



associated with the observable a{x) captures the dynamics of both the flow and the 
Poincare map. The dynamical zeta function for the flow is obtained as l/^{s,fi) = 
1 /^(l, s,fi), and the dynamical zeta function for the discrete time Poincare map is 

1/^(Z,^)- 1/^(Z,0,;S). 

Our basic assumption will be probabilistic. We assume that the chaotic in- 
terludes render the consecutive return (or recurrence) times T(xj), r(x,+i) and ob- 
servables a(xi), a(xi+i) effectively uncoiTclated. Consider the quantity gPMxo,n)-sT(xo,n) 
averaged over the surface of section P. With the above probabilistic assumption 
the large n behavior is 



^^A(xo.n)-sT(xo.n)y^ ~ |J e^"(^)-^V(;c)t/;cJ , 



where p(x) is the invariant density of the Poincare map. This type of behavior is 
equivalent to there being only one zero zo(^,l3) = j e^"^^^''^'^^^^p{x)dx of 1 /^(z, s,^) 
in the z-/3 plane. In the language of Ruelle-PoUicott resonances this means that 
there is an infinite gap to the first resonance. This in turn implies that 1 /^(z, s,/3) 
may be written as remark 17.1 



l/az,s,/S) = z- (f'''''-'^'''p{x)dx , (24.44) 



where we have neglected a possible analytic and non-zero prefactor. The dynam- 
ical zeta function of the flow is now 



\l^{s,f3) = \l^{\,s,l5) = \- r ^''^'''^p{x)e-'^^'''^dx . (24.45) 

Jp 

Normally, the best one can hope for is a finite gap to the leading resonance of 
the Poincare map. with the above dynamical zeta function only approximatively 
valid. As it is derived from an approximation due to Baladi, Eckmann, and Ruelle, 
we shall refer to it as the BER zeta function 1/^ber('^'/3) iri what follows. 

A central role is played by the probability distribution of return times 

iA(t) - r 6{t- T{x))p(x)dx (24.46) 
Jp 



exercise 25.6 
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The BER zeta function at y6 = is then given in terms of the Laplace transform of 
this distribution 

1/^BEr(^) = 1- Hr)e'''dT. 
Jo 

exercise 24.5 

Example 24.1 Return times for the Bernoulli map. For the Bernoulli shift map 
(23.6) 

X i-» f{x) - 2x mod 1 , 

one easily derives the distribution of return times 
1 

The BER zeta function becomes (by the discrete Laplace transform (18.9)) 

CO CO ^ 

«- 1 n- 1 

= = r'W/d-z/Ao). (24.47) 

1 - z/2 

Thanks to the uniformity of the piecewise linear map measure (16.13) the "approximate^ 
zeta function is in this case the exact dynamical zeta function, with the periodic point 
pruned. 



Example 24.2 Return times for the model of sect. 24.2.1. For the toy model of 
sect. 24.2.1 one gets ip\ - \M\\, and tfr„ = |A1„I(1 - b)/{l - a), for n > 2, leading to a 
BER zeta function 

DO 

n=2 

which again coincides with the exact result, (24. 11). 

It may seem surprising that the BER approximation produces exact results in 
the two examples above. The reason for this peculiarity is that both these systems 
are piecewise linear and have complete Markov partitions. As long as the map 
is piecewise linear and complete, and the probabilistic approximation is exactly 
fulfilled, the cycle expansion curvature terms vanish. The BER zeta function and 
the fundamental part of a cycle expansion discussed in sect. 20.1.1 are indeed 
intricately related, but not identical in general. In particular-, note that the BER zeta 
function obeys the flow conservation sum rule (20. 17) by construction, whereas 
the fundamental part of a cycle expansion as a rule does not. 
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Resume 

The presence of marginally stable fixed points and cycles changes the analytic 
structure of dynamical zeta functions and the rules for constructing cycle ex- 
pansions. The marginal orbits have to be omitted, and the cycle expansions 
now need to include families of infinitely many longer and longer unstable or- 
bits which accumulate toward the marginally stable cycles. Correlations for such 
non-hyperbolic systems may decay algebraically with the decay rates controlled 
by the branch cuts of dynamical zeta functions. Compared to pure hyperbolic 
systems, the physical consequences are drastic: exponential decays are replaced 
by slow power-law decays, and transport properties, such as the diffusion may 
become anomalous. 



Commentary 

Remark 24.1 What about the evolution operator formalism? Tlie main virtue of 
evolution operators was their semigroup property (17.27). Tfiis was natural for hyper- 
bofic systems wliere instabifities grow exponentiaffy, and evofution operators capture ttiis 
behavior due to their muftipficative nature. Whether the evofution operator formafism is a 
good way to capture the sfow, power f aw instabifities of intermittent dynamics is fess cfear. 
The approach tatcen here feads us to a formuf ation in terms of dynamical zeta functions 
rather than spectraf determinants, circumventing evofution operators aftogether. ft is not 
fcnown if the spectraf determinants formufation woufd yiefd any benefits when appfied to 
intermittent chaos. Some resufts on spectraf determinants and intermittency can be found 
in [24.2]. A usefuf mathematicaf technique to deaf with isofated marginaffy stabfe fixed 
point is that of inducing, that is, repf acing the intermittent map by a compfetefy hyperbofic 
map with infinite afphabet and redefining the discrete time; we have used this method im- 
pficitfy by changing from a finite to an infinite afphabet. We refer to refs. [24.3, 24.20] 
for detaifed discussions of this technique, as weff as appfications to 1-dimensionaf maps. 

Remark 24.2 Intermittency. Intermittency was discovered by Manneviffe and Pomeau [24. 1] 
in their study of the Lorenz system. They demonstrated that in neighborhood of param- 
eter vafue = 166.07 the mean duration of the periodic motion scafes as (r - r^)''^. 
In ref. [24.5] they explained this phenomenon in terms of a 1-dimensionaf map (such as 
(24. 1)) near tangent bifurcation, and classified possibfe types ofintermittency. 

Piecewise finear modefs Ufce the one considered here have been studied by Gaspard 
and Wang [24.6]. The escape probfem has here been treated foffowing ref. [24.7], resum- 
mations foffowing ref. [24.8]. The proof of the bound (24.28) is given in P. Dahfqvist's 
notes, see ChaosBook.org/extras/PDahlqvistEscape.pdf. 

Farey map (24.42) has been studied widefy in the context of intermittent dynamics, 
for exampfe in refs. [24.16, 24.17, 20.3, 24.18, 24.19, 20.16, 24.2]. The Fredhofm deter- 
minant and the dynamicaf zeta functions for the Farey map (24.42) and the refated Gauss 
shift map (16.47) have been studied by Mayer [24.16]. He relates the continued fraction 
transformation to the Riemann zeta function, and constructs a Hifbert space on which the 
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evolution operator is self-adjoint, and its eigenvalues are exponentially spaced, just as for 
the dynamical zeta functions [24.24] for "Axiom A" hyperbolic systems. 



Remark 24.3 Tauberian theorems. In this chapter we used Tauberian theorems for 
power series and Laplace transforms: Feller's monograph [24.9] is a highly recommended 
introduction to these methods. 



Remark 24.4 Probabilistic metliods, BER zeta functions. ProbabiHstic description 
of intermittent chaos was introduced by Geisal and Thomae [24.10]. The BER approx- 
imation studied here is inspired by Baladi, Eckmann and Ruelle [24.14], with further 
developments in refs. [24.13, 24.15]. 
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Exercises 



24.1. Integral representation of Jonquiere functions. 

Check the integral representation 



24.6. 



24.4. 



24.5. 



J{z, a) 



r(a) 



f 

Jo 



for a > .(24.48) 



Note how the denominator is connected to Bose- 
Einstein distribution. Compute J{x + ie) - J(x - ie) for 
a real x> 1 . 

24.2. Power law correction to a power law. Expand 
(24.21) further and derive the leading power law correc- 
tion to (24.24). 

24.3. Power-law fall off. In cycle expansions the stabilities 
of orbits do not always behave in a geometric fashion. 
Consider the map / 




This map behaves as / — > x as — > 0. Define a sym- 
bolic dynamics for this map by assigning to the points 
that land on the interval [0, 1 /2) and 1 to the points that 
land on (1/2, 1]. Show that the stability of orbits that 
spend a long time on the side goes as In particular, 
show that 

Power law fall-off of Floquet multipliers in the sta- 
dium billiard**. From the cycle expansions point 
of view, the most important consequence of the shear in 
J" for long sequences of rotation bounces in (8.13) 
is that the A„ grows only as a power law in number of 
bounces; 



A„ oc «,, 



(24.49) 



Check. 



Probabilistic zeta function for maps. Derive the 
probabilistic zeta function for a map with recurrence dis- 
tribution 41 „. 



24.7. 



Accelerated diffusion. Consider a map h, such that 
h - f, but now running branches are turner into stand- 
ing branches and vice versa, so that 1,2,3,4 are stand- 
ing while leads to both positive and negative jumps. 
Build the corresponding dynamical zeta function and 
show that 



t 

tint 

fVlnf 
,2 



for a > 2 
for a - 2 
for ae (1,2) 
for a - I 
for a € (0, 1) 



Anomalous diffusion (hyperbolic maps). Anoma- 
lous diffusive properties are associated to deviations 
from linearity of the variance of the phase variable we 
are looking at: this means the diffusion constant (17.13) 
either vanishes or diverges. We briefly illustrate in this 
exercise how the local local properties of a map are cru- 
cial to account for anomalous behavior even for hyper- 
bolic systems. 

Consider a class of piecewise linear maps, relevant to 
the problem of the onset of difiiision, defined by 



Ax 

a — Af^ylx - 
1 - A'(x-x+) 
l-a 
1 H- A(x- 1) 



Af^ylX - X 



for 


X e 


0, 


4] 


for 


X e 




,x- 


for 


X e 


x+ 

2 


,x- 


for 


X € 




, X2 


for 


X 6 


X2 


■■1 



where A = (1/3 - e'^'')-'. A' = (1/3 - 2e^'^), A,,y = 
e'-i^^fl = l+e,x+ = l/3,x| =x+-6i/'',x+ ^x+H-e'^^ 
and the usual symmetry properties (25. 1 1) are satisfied. 

Thus this class of maps is characterized by two escap- 
ing windows (through which the diff'usion process may 
take place) of size 2e'^'': the exponent y mimicks the or- 
der of the maximum for a continuous map, while piece- 
wise linearity, besides making curvatures vanish and 
leading to finite cycle expansions, prevents the appear- 
ance of stable cycles. The symbolic dynamics is eas- 
ily described once we consider a sequence of param- 
eter values {e,„), where e,„ - A"'"'^'': we then par- 
tition the unit interval though the sequence of points 
0, x|, x^, Xj, Xp x", xj, 1 and label the corresponding 
sub-intervals 1, Sa, Sb,2,db,da,3: symbolic dynamics is 
described by an unrestricted grammar over the following 
set of symbols 



{l,2,3,i#- l'',^f#- 3* 



# — a,b i,k — m,m + \,m 
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This leads to the following dynamical zeta function: a matter of fact, from (24.51) we get the asymptotic be 

2^ ^ jhavior D ~ e'^'', which shows how the onset of diffusioi 

^0 '(z, a) - ^ ~ ~^ ~ ^ ~ ^ cosh(a)e^^''"' ~ x) governed by the order of the map at its maximum. 



from which, by (25.8) we get 
D = 



Remark 24.5 Onset of diffusion for continuous maps. 

2e,J/'' 'A"'"(1 - 1/A)"' The zoology of behavior for continuous maps at the 

I _ 2_ 1 4f'/''"' / m+ 1 ^ 1 \^^'^'^^^nset of diffusion is described in refs. [25.15, 25.16, 

^ '" VA"'(i-i/A) A"'+ (i-i/A)-/ 24.25]: our treatment for piecewise linear maps was in- 

The main interest in this expression is that it allows ex- troduced in ref. [24.26]. 

ploring how D vanishes in the e (m i-> oo) limit: as 
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